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Abstract: We propose a method to assess the quality of geometrical scaling in Deep 
Inelastic Scattering and apply it to the combined HERA data on j*p cross-section. Using 
two different approaches based on Bjorken x binning and binning in j*p scattering energy 
W, we show that geometrical scaling in variable r ~ Q 2 x x works well up to Bjorken x's 
between 0.02 and 0.08. The corresponding value of exponent A is 0.32 - 0.33. 
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1 Introduction 

Geometrical scaling (GS) introduced in Ref. [1] in the context of low x Deep Inelastic 
Scattering (DIS) is by now well established phenomenon attributed to the existence of 
an intermediate energy scale, called saturation scale Q s (x). Saturation scale is defined 
as the border line between dense and dilute gluonic systems within a proton (for review 
see e.g. Refs. [2, 3]). GS has been also observed in DIS on nucleus [4] and in diffraction 
[5]. Recently it has been shown that GS is also exhibited by the pr spectra at the LHC 
[6-8]. Geometrical scaling has a natural explanation within the theory of saturation and 
the Color Glass Condensate [9-12]. 

Geometrical scaling takes place if some observable that in principle depends on two 
independent kinematical variables like Q 2 (or p 2 -,) and W (i.e. ^*p c.m.s. energy) depends 
only on a specific combination of them, namely on 



Q 2 



called scaling variable. Here 



Q 2 



(l.i) 



(1.2) 



Q2 + W 2_ M 2 

is Bjorken x variable, M p stands for the proton mass, and the saturation scale Q s takes 
the following form [1, 13] 

Ql(x) = Ql (±) . (1.3) 



Here Qo and xq are free parameters which can be fitted to the data within some specific 
model of DIS, and exponent A is a dynamical quantity of the order of A ~ 0.3. Phys- 
ical observable that exhibits GS in the case of DIS is 7*p cross-section <jy* p (x, Q 2 ) = 
47r 2 a em F 2 (x,Q 2 )/Q 2 . 

Although many authors, following the original paper by Stasto, Golec-Biernat and 
Kwecinski [1], have shown that GS is seen in DIS data, until now there was no quantitative, 
model independent analysis of its applicability domain. Different forms of scaling variable 
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have been tested in a series of papers [14-16] where the so called Quality Factor (QF) has 
been defined and used as a tool to assess the quality of geometrical scaling. These authors, 
however, constrained their analysis only to the domain of small Bjorken x's, x < 0.01 and 
concentrated on testing the different forms of scaling variable r. The QF is a new tool 
for which the quantitative statements, like the one concerning the value of exponent A for 
example, are not based on the standard chi-square analysis. 

In contrast, in this paper we propose a new numerical criterion for GS which serves as 
a tool to extract exponent A by standard \ 2 minimization. We consider cr^*p{xi, Q 2 ) for 
different fixed functions of Q 2 . Geometrical scaling hypothesis means that 

a rp ( Xi ,Q 2 ) = -Lf(t) (1.4) 

where F(r) is a universal dimensionless function of r. Therefore if cross-sections a~* p {xi, Q 2 ) 
for different Xi are evaluated not in terms of Q 2 but in terms of r, they should fall on one 
universal curve (see Figs. 11 and 12 in Sect. 4). This in turn means that if we calculate 
ratio of cross-sections for different Bjorken Xi's each expressed in terms of r, we should get 
unity independently of r. This allows to determine power A by minimizing deviations of 
these ratios from unity. We can form as many ratios as there are different pairs of Xj's in 
overlapping regions of Q 2 (or more precisely of r). Needless to say that the best values of 
A's extracted from different ratios should coincide within errors. 

In an ideal case one would choose the lowest possible x as the reference x re f to calculate 
the ratios of the cross-sections. Unfortunately there is a strong correlation between Bjorken 
x's and values of Q 2 measured by HERA [17], and therefore there is no single value of X{ 
which covers all available values of Q 2 . It turns out that the coverage in space of W (rather 
than x) and Q 2 is much broader. Therefore in what follows we also study the quality of 
GS for a^piWi, Q 2 ) in bins of W although this requires "rebinning" of the data which are 
provided by HERA experiments only in bins of (x,Q 2 ). There is, however, an advantage 
of such a procedure, as it is almost identical to the analysis applied to the pr spectra in 
pp collisions at the LHC [6, 8]. 

Since our analysis is sensitive only to the variations of scaling variable r with x and not 
to the absolute value of r, we choose in the following (unless specified otherwise) Qo = 1 
GeV/c and xq = 1, i.e.: 

T = Q 2 x x . (1.5) 

The absolute value of the saturation scale can be inferred only from some explicit model 
of DIS at low x. For the purpose of the subsequent analysis we define "y*p cross-section as 

, n 2, F 2 (x,Q 2 ) 
a rp (x,Q ) = — 2 (1.6) 

since all proportionality constants cancel out in the ratios of the cross-sections that are the 
main tool used to look for geometrical scaling in this paper. The corresponding experi- 
mental error is therefore given by 

A<^(^) 2 + (§AQf. d.7) 
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HERA kinematic range 
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Figure 1. HERA kinematic range in (Q 2 ,x) plane. It is shown which Bjorken x's (brown full 
squares) are chosen for given x re f (red stars) to construct ratios of ~f*p cross-sections in scaling 
variable r. 

For AF2 we use the measured errors of the structure function, for Q 2 we calculate the error 
from the widths of the Q 2 bins (for details see Ref. [18]). 

Our findings can be shortly summarized as follows. Geometrical scaling in variable 
(1.5) works well up to Bjorken x's between 0.02 and 0.08. The corresponding value of 
exponent A is 0.32 - 0.33. 

The paper is essentially divided into two parts. In Sect. 2 we discuss GS in terms of 
Bjorken variable x, whereas in Sect. 3 in terms of 7*2? scattering energy W, in both cases 
for e + p scattering. In Sect. 4 we summarize our results and discuss other experiments and 
also HERA data for e~p. We also present possible applications of the method proposed in 
this paper to test other possible forms of scaling variable r and to investigate scaling of 
the charm cross-section. 

2 Bjorken x binning 

Combined study of ZEUS and HI [17] provides 432 e + p data points for DIS structure 
function F2 in terms of Bjorken x and Q 2 . Here we shall use 348 points with Bjorken x's 
that have at least 2 points in Q 2 [18]. This choice corresponds to 5.52 x 10 -6 < x < 0.65. 
Since there is a strong correlation between x and Q 2 as far as HERA kinematic range 
is concerned (see Fig. 1), the analysis outlined in Sect.l requires some refinement. We 
apply here the following procedure. First we choose some x re f and consider all Bjorken 
Xj's smaller than x re f that have at least two overlapping points in Q 2 (or more precisely 
in scaling variable r), as depicted in Fig. 1 (note, however, that when we make the same 
plot in terms of variable r rather than Q 2 the overlaps change, usually the number of 
overlapping points gets smaller). We require that x» < x re f because GS holds for small 
values of Bjorken x and therefore by increasing x re f we will be able to see violations of GS. 
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For Xi < x re f the ratios 

R X i,x ief (^, Tk) = U ] P Xl,T ( Xt, ®k' X )) with Tk = T ^ q2. ^ = r ( Xref) Q2 A) (2.1) 



<J-y* 



(x rei ,T(x rei ,Q 2 kre{ ;\)) 



are greater than 1 for A = 0. By increasing A one shifts R Xi ,x iaf (A; Tk) towards unity for all 
values of Tk- 

Note that t{xi,Q\\X = 0) = Q\ and ratios (2.1) are essentially ratios of 7*p cross- 
sections at overlapping values of Q\ (as in Fig. 1). For A / points of the same Q 2 
but different x's correspond generally to different r's. In order to calculate ratios (2.1) at 
points Tk corresponding to one fixed Xi, one has to interpolate the reference cross-section 
p-y*p{x ie f, T(x re f, Q 2 ; A)) to Q| ref such that r(x re f, Q\ ref ; A) = Tfc. Since 7*p cross-sections 
to a good accuracy lie on straight lines as functions of logQ 2 , we apply in the following 
the linear interpolation in logQ 2 . Interpolation errors are neglected. 

It is clear than not all available Bjorken x's can be used as x re f. If one chooses x Te { 
too small then there are no x < x re f which have at least two points in r (or Q 2 ) within 
cr 7 * p (r(x re f, Q 2 )) domain. It turns out that we can use only x re f > 3.2 • 1CT 5 . 

In order to find optimal exponent A that minimizes deviations of ratios (2.1) from 
unity we form the chi-square measure 



> (\) = 1 (^,x rcf (A; r fc ) - l) 2 



where the sum over k extends over all points of given Xi that have overlap with x Te {. As 
already explained above, for A = this is essentially the sum over all overlapping values of 
Q 2 , for A / we choose measured points of o" 7 * p (r(xj, Q\; A)) and interpolate the reference 
cross-section in Q 2 to the point Q 2 rof such that r(x re f, Q\ ref ; A) = Tk- 
Finally, the errors entering formula (2.2) are calculated using 



ARxi,x rcf {X]Tk) 



I Aa rp (xi,T(xi,Q 2 k )) \ 

I (J^*p(x Te {, T(x re f, Q 2 re f)) I 

where /S.a^* p {r(x, Q 2 )) are experimental errors (or interpolated experimental errors) of *y*p 
cross-sections (1.7). 

By minimizing X x ,x re fW °f -^l- (2-2) with respect to A one can compute the "best" 
value A m i n (x, x re f) and its error requiring that {N XuXtc{ — l)x Xi ,x Tcf (^rnm) changes by 1 around 
the minimum. In the region where GS is satisfied A m i n should not depend neither on x nor 
on x re f and Xx,z ref (Amin) should be small. In Fig. 2 we plot ratios R x ,x Ici (X;Tk) as functions 
of Tk for A = and A = A m i n for some typical values of x and x re f. In Fig. 2. a we plot R 
for both x and x re f small and quite distant, i.e. x re f/x = 10. Here we expect GS to be 
satisfied. We see that for A = ratio R grows with Tk = Q\ and is of the order 1.5 -r- 2. 
By minimizing Xx,x ie{ (^) we obtain A m i n = 0.306, which agrees with the expectations. In 
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Figure 2. Ratios Ri(\) for A = and A = A m j n for some typical choices of x and x Ie {- 



Fig. 2.b, both x and x re f are small, however they are not very distant: x ve {/x = 1.6. We see 
that already for A = ratio R XiXie{ (0; r^) is close to 1; we can improve the value of x 2 by 
increasing the value of A to A m i n = 0.308, in fair agreement with the previous case, however 
one should note that now X 2 W is rather flat. Another example is shown in Fig. 2.c where 
%ref = 0.25 is relatively large and x still quite small, such that x Te f/x = 5. Here again it is 
possible to make i? XiXret (A; T&) close to unity, however at the expense of rather large value 
of A m i n = 0.44. This is a clear sign of violation of the universality of exponent A, which we 
consider to be a signature of violation of geometrical scaling. Similarly in Fig. 2.d we show 
Rx,x re{ (A; Tfc) for both x and x re f large and not so close: x m {/x = 2.22. Here A m ; n = 0.965, 
which is outrageously large. 

The qualitative measure of geometrical scaling is given by the independence of A m i n 
on Bjorken x and by the value of X x ,x rei (^min) in the minimum. In Fig. 3 w show contour 
plots of A m i n (x, x re f) and X x ,x Tct (^min) in the (x, x re f) plane. By construction both A m i n and 
X 2 (^min) are defined only above the line rr re f = x. We see From Fig. 3. a that the stability 
corner of A m i n (i.e. the region where variations of A m j n are small ) extends up to x re f < 0.2. 
In the most of this region 0.3 < A m i n < 0.4. Similar pattern can be seen in Fig. 3.b where 
X x , Xiei (Amin) starts rising above x re f ~ 0.35 with some fluctuation below. From Fig. 3 we 
can conclude that geometrical scaling holds up to Bjorken x's of the order of 10 _1 which 
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Figure 3. Contour plots of X m - m (x,x IC [) (a) and xi x, ot (Amin) (b) in the x,x IC [ plane. Level lines 
are spaced by 0.1 for A and by 0.15 for \ 2 - To set the scale two sample level lines are labeled: 0.4 
for A and 0.15 for x 2 . 



is well above the original expectations. This is nicely illustrated by the three dimensional 
plot of A m i n (x, x Te {) in Fig. 4 where the scale of the axes x and x re f is logarithmic. We see 
rather flat distribution of A m i n up to x re f ~ 10" 1 and rapid increase for larger Bjorken x's. 

We can now look for possible violation of GS in more quantitative way. In order to 
eliminate the dependence of A m i n (x, x re f) on the value of x, we introduce averages over x 
(denoted in the following by (...)) minimizing the following chi-square function: 



2 



~2 ^\\^ 1 {^mm\X, X Ie f) (A)) . . 

**,.,(«) = jt— T _£ AKUx ^, v (2-4) 

which gives the "best" value of A denoted as (A m i n (x re f)). The sum in (2.4) extends over 
all x's such that A m ; n (x, x re f) exists (see Fig. 1). N Xiel is the number of terms in (2.4). 
The results are plotted in Fig. 5 as full (black) squares. Errors A (A m i n ) are calculated 
from the requirement that (N Xiei — l)x^ ref (A m in) changes by 1. For some points one can 
see very large errors; this is due to the small number of possible x's for this particular x re f. 
On average, however, one can approximate (A m i n (x re f)) with a constant value of 0.32 up 
to x TC { ~ 0.1 where the rise of (A m i n (x re f)) can be interpreted as violation of geometrical 
scaling. 

Looking at Fig. 5 one may have an impression that on the average (A m i n (x re f)) is 
slightly rising with x. This small rise might be, however, attributed to the dependence of A 
on Q 2 through the kinematical correlation between Bjorken x's and Q 2 in the kinematical 
HERA range (see. Fig. 1). Q 2 dependence of A has been measured in DIS for very small 



Figure 4. Three dimensional plot of X m i n (x,x rc f). 
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Figure 5. Left: averaged values (A m i n (a; re f)) (black squares) and ((A m j n (a; cu t))) (red circles); right: 
Xx { (black squares) and Xz cut ( re d circles) as functions of x re i and x cu t, respectively. 



x's [19] and has been discussed in the context of the pt spectra at the LHC [7]. It is 
theoretically motivated by the corrections due the DGLAP evolution [20, 21]. We shall 
come back to this possibility in Sect. 4. 

To quantify further the hypothesis of geometrical scaling with constant A we form yet 
another chi-square function 



-2 f//\\\\ 1 (^min(x, X Te f) — ((A))) 

Cx cut {\\A))) - k — r 2- 

A " x ^ x rcf <a 

which we minimize to obtain ((A m i n (rc cu t))). 



Nx <«* ~ 1 x ttx *x<?, AA min (x,a; ref ) 2 
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Figure 6. HERA kinematic range in (Q 2 , x) plane with several lines of constant W. 

The idea behind equation (2.5) is to see how well one can fit (A m i n (x re f)) with a constant 
A up to x re f = x cut . Were there any strong violations of GS above some xq, one should 
see a rise of ((A m i n (x cut ))) once x cut becomes larger than xq. As can be seen from Fig. 5, 
where ((A m i n (x cu t))) and the corresponding Xx cut are plotted as full (red) circles, no drastic 
change in ((A m i n (x cut ))) can be seen, although the slow rise is clearly seen above x cut ^ 0.1. 
On the contrary Xx cut r i ses steeply above x cut ^ 0.1. Since for large x cut the number of 
points N Xcut is of the order of a few tens, we can safely assume that violation of GS takes 
place when Xx cut > 0.3 which happens precisely around x cu t — 0.1. We find therefore that 
the best value of a constant A corresponds to x cu t = 0.08 which we denote as 

A Bj = 0.319 ± 0.005 (2.6) 

where subscript "Bj" stands for Bjorken x binning. 

Let us finish by a remark on comparison of ((A m i n (x cu t))) with a more "differential" 
quantity (A m i n (x re f)) depicted in Fig. 5. Since ((A m i n (x cut ))) is in a sense an average 
of all (A m i n (x re f)) for x re f < x C ut, the rapid increase of (A m i n (x re f)) above x re f ~ 0.1 is 
smoothed out due to a long constant tail of (A m i n (x re f)) for smaller x re f's. This happens, 
however, at the expense of the increase of the corresponding x 2 - Below 0.1 we have that 
(A m i n (x)) ~ ((A m j n (a;))) which further confirms validity of geometrical scaling in this region. 

3 Energy binning 

Here we are going to repeat analysis of Sect. 2 for the combined HERA data [17], but in 
bins of (W, Q 2 ) rather than (x, Q 2 ). The reason is that fixed energy points span over much 
wider common Q 2 range (see Fig. 6) and therefore we do not have to to choose different 
reference points and average over them. Moreover, analysis in terms of W follows exactly 
the method used to analyze pt spectra at the LHC [6, 8]. The disadvantage consists in the 
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Table 1. Energy bins and energies assigned to them. Number of points in different bins is also 
displayed (these are values for e + p data). 



fact that HERA points are not organized in the energy bins, therefore we have to "rebin" 
them using formula (1.2). This leads to the loss of statistics. Indeed, instead of 432 points 
(see Sect. 2) we shall use here 347 e + p data points. Furthermore some ambiguities arise, 
such as the choice of the bin size, the need to recalculate Bjorken x and others discussed in 
detail in Ref. [18]. In what follows we shall use logarithmic binning with step 1.3 — every 
consecutive W is 1.3 times larger than preceding one (in Ref. [18] different binnings have 
been also considered without any major change of the results). We define value of energy 
W as the mean of two limiting values W^ in and W4ax between which it lies (see Table 1). 

In what follows we use W re f = 206 GeV because it gives the widest range of r values 
and is one of the biggest energies that we have at our disposal (GS is expected to be 
present at large energies), and it has the largest number of points in Q 2 . In Ref. [18] 
different choices of W re f have been also analyzed with essentially the same conclusions. 

Similarly to the case of x -binning we shall construct ratios of cross-sections at 
different Ws as functions of r. In general r values for W re f and W are different so we 
need to interpolate the reference cross-section to the value of Q 2 c{ which corresponds to the 
r value needed to calculate the ratio. Here we follow closely the method of Sect. 2 using 
linear interpolation in log Q 2 and neglecting interpolation errors. 

Now for every point with energy Wi (note that W% < Wref) we define a ratio (k labels 
points with energy W% but of different Q 2 ): 

R Wi,w iel ( X , T k) ■= 7 (w n vy, W1 th r k = T (Wi, Q k ; A) = r(W ref , Q k f ; A). 

(3-1) 

Uncertainty of ratio Rw it w iei (^] T k) is given by: 

/ Aa rp (W vci , r(W rei , Q 2 ref ))\ 2 



A-RWi.WrefCAjTfc) 



P WMWi,Ql)) 



+ { v rp m,TW,Q 2 )) ^WMWiM)] . (3.2) 
Since we have fixed W m f = 206 GeV, in what follows we shall omit subscript Wr e f. 
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Figure 7. Ratios for A = and A = A m ; n = 0.367 both with x cut = 0.25. 



In Fig. 7 we show as an example R72 plotted for A = and A = 0.367. We see that 
R72 decreases several times when we use scaling variable rather than Q 2 . This is generic 
feature which we employ to look for geometrical scaling. 

Our aim is to find such A for given energy W% < W re f that deviations Ry/ t (A; T k ) — 1 
are minimal. Taking into account uncertainties ARw we define the chi-square function: 



XWi ( x cut', A) 



1 



E 



(R Wi (\]T k ) 



I) 



1 ^ AR w .(\;r k )2 

keWr, x<x cnt 



(3.3) 



where k £ Wf, x < x cn t means that we sum over points corresponding to given energy Wi 
and values of x that are not larger than x cn t- The reason to introduce a cut-off on Bjorken 
x is to look for violations of GS once we get into the region of large x's. 

We will search A m i n (Wi, x cut ) which minimizes Xw- ^ or §i ven W% and x cut . Uncertainty 
of A m ; n is estimated by requiring that (N\y z ,x cut — l)Xw- ( x cut] A) changes by 1 when A is 
varied around the minimum. The results for 12 different energy bins are plotted in Fig. 8. 
One can see that in each Wi bin A m i n (x cu t) can be approximated by a constant, although for 
lower energies a slight increase for x cu t > 0.1 is present. Large error bars for higher energies 
are due to the fact that %w- S e ^ s flatter once Wi is close to W re f. Since at low energies there 
are no points of low x < 10 -3 , in what follows we shall restrict the energy range to Wi > 25 
GeV. We shall also skip the highest energy bin of W = 268 GeV since it contains only 7 
points and the corresponding A m i n has very large errors (see Fig. 8) exhibiting different 
behavior with x cu t than for other energies, which suggests some systematic binning effect 
that we do not have under control. Therefore the final number of data points used in the 
analysis in this Section is further reduced to 318. In Fig. 9 we plot the corresponding values 
of Xw - O ne can see that up to x cut ^ 0.1 the X%/^ s are below 0.2, and the rapid growth 
is seen for x cut > 0.1. Before we can conclude that GS is well satisfied for Bjorken x's up 
to 0.1, we have to check whether the corresponding X m i n (Wi) are energy independent. For 
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Figure 8. Exponent A m j n as a function of a; cut for all energies W ^ W rc f = 206 GeV. 



very small x cu t's in turn we have small numbers of points in (3.3) and the calculated values 
of X\v are statistically not reliable. 

Finally, in order to get rid of the energy dependence we define yet another chi-square 
function 

~1i \\ 1 (Amin (Wi,X cut) ~ A) 2 

X (zcut; A) = — > — — — ^— (3.4) 



N 



^cut 



AA min (Wi 



where the sum goes over eight energies Wi = 25, 33, . . . , 159 GeV. By minimizing (3.4) 
with respect to A we obtain the best value of A m i n (Wi, x cu t) averaged over all energies 
which we denote by A ave (x cu t). The error of A ave is calculated again by demanding that 
(-^scut — l)x 2 ( x cut] A) changes by 1 when A is varied around A ave . Both Aave(^cut) and 
the corresponding x 2 are plotted as functions of x cu t in Fig. 10. We see that A aV e(^cut) is 
rather flat; large errors for small x cu t are due to the small number of points with small x's. 
Looking at A ave (x cu t) one does not see any dramatic change for x cut — > 1. On the other 
hand the value of x 2 starts rapidly growing at x cu t — 0.02. We can therefore conclude that 
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Figure 10. Exponent A m i n averaged over energies (denoted as A avo ) as a function of x cut and the 
corresponding % 2 . 



GS works well up to Bjorken x's of the order of 0.02 with 

A en = A ave (x cut = 0.02) = 0.328 ± 0.011. (3.5) 
where subscript "en" stands for energy binning. 

4 Summary and outlook 

In the present paper we have performed quantitative, model independent analysis of the 
accuracy and the applicability domain of geometrical scaling in deep inelastic e + p scatter- 
ing. To this end we have chosen the most recent compilation of the HERA data based on 
common analysis of ZEUS and HI experiments [17]. We have tested the standard form of 
geometrical scaling [1, 13] given by the form of the scaling variable r (1.5), assuming the 
constant value of the exponent A. In order to quantify the quality of GS we have proposed 
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two different approaches. In the first approach we have computed ratios of cross-sections 
F2/Q 2 for two different Bjorken re's, denoted as x and rc re f, as functions of the scaling 
variable r; such ratios should be equal to unity if GS is present. This allowed us to extract 
the best value of the exponent A, which we called A m i n (a;, x re f). In the region where GS is 
present A m i n (x, x re f) should not depend on x and x re f, and moreover the corresponding x 2 
should be small. In the second approach we have repeated essentially the same calculation 
but for 7*p cross-sections treated as functions of the scattering energy W rather than x. 
This procedure is analogous to the one applied to the pr spectra at the LHC. It requires, 
however, "rebinning" of the data, which reduces the statistics and introduces uncontrol- 
lable errors, for example the uncertainty of x, and also in principle some uncontrollable 
dependence on the choice of W bins. Nevertheless, since "f*p cross-sections for all available 
energies overlap in Q 2 (or r) nearly over whole kinematical range (which is not the case 
for the x binning, see Figs. 1 and 6), the analysis is much more straightforward than in 
the case of the Bjorken x binning. 

Despite the differences between the Bjorken x binning and the energy binning the 
results in both cases are consistent and to some extent surprising. By averaging A m i n (x, x re {) 
or X m in(W, W re f) over the allowed regions of x or W respectively, and by inspecting the 
corresponding x 2 values, as described in detail in Sects. 2 and 3, we have established that 
geometrical scaling holds up to relatively large x's of the order of 0.08 for the Bjorken x 
binning and 0.02 for the energy binning. The resulting exponents A given by Eqs. (2.6) and 
(3.5) agree within errors and are consistent with previous estimates and model calculations 

A = 0.32 -h0.33. (4.1) 

These results are suggestively illustrated in Figs. 11 and 12 where we plot 7*p cross-sections 
for both types of binning used in this paper as functions of Q 2 and the scaling variable r 
(with Qq = 1 GeV/c and xo = 10~ 3 in Eq.(1.3)). One can see that in both cases the cross- 
sections fall on one universal curve - a clear sign of geometrical scaling. Some deviations 
of GS can be seen in the lower right end of the plots. These points, however, correspond 
to large x's (due to the kinematical correlation of the HERA phase space, see Figs. 1,6), 
and therefore should not exhibit geometrical scaling. 

The analysis performed here is essentially identical to the one presented in in Ref. [18] 
with one important difference. Here, when interpolating cross-section in Q 2 we used 
linear interpolation in logQ 2 , rather than in Q 2 . The latter interpolation introduces large 
errors which were neglected in [18]. This resulted in somewhat different values of the 
exponent A. Qualitative picture, however, remains unchanged. 

One of the questions addressed in [18] was possible dependence of A on Q 2 . Apply- 
ing the same method of dividing j*p cross-sections of different Bjorken x's (or scattering 
energies W) we have been looking for the best exponent A which, in this case however, 
depended on Q 2 . Instead of having one parameter A we had therefore a number of pa- 
rameters Afc where k runs over all Q 2 bins. By minimizing x analogous to (2.2) or (3.3) 
we did not observe any visible dependence of A on Q 2 . This is most probably due to the 
fact that the corresponding x 2 functions are very flat and one can easily fall into some 
accidental minimum. We therefore conclude that with the present experimental accuracy 
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Figure 11. Geometrical scaling for the Bjorken binning. Left: j*p cross-sections F2/Q 2 as func- 
tions of Q 2 for fixed x. Different points correspond to different Bjorken x's. Right: the same but 
in function of scaling variable r for A = 0.319. Points in the right end of the plot correspond to 
large x's (due to kinematical correlation of the HERA phase space), and therefore show explicitly 
violation of geometrical scaling. 
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Figure 12. Geometrical scaling for the energy binning. Left: 7*p cross-sections F2/Q 2 as functions 
of Q 2 for fixed W's. Right: the same but in function of scaling variable r for A = 0.328. Points 
in the right end of the plot correspond to large x's and mostly small energies (due to kinematical 
correlation of the HERA phase space), and therefore show explicitly violation of geometrical scaling. 



we have been unable to find Q 2 dependence of A in a model independent way. It seems 
that in this case one has to assume some Ansatz for \{Q 2 ) depending on a few variational 
parameters which can be found using the methods described in this paper. Universality of 
these parameters would then constitute a criterion for Q 2 dependence of A. 

We have also analyzed separately j*p cross-section for e~p scattering. The e~p data 
cover narrower range of Q 2 values and therefore the smallest x mm = 0.002 is rather large 
as far as GS is concerned. Nevertheless the results for e~p and e + p scattering restricted to 
the same kinematical region are very similar except for some small systematic shift of the 
exponent A towards higher values in the e~p case [18]. 

Natural question concerning other deep inelastic scattering experiments arises. Here 
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we have \ip scattering experiments: EMC [22] with the smallest Bjorken x m { Q = 0.0175, 
BCMDS [23] (x min = 0.07), NMC [24] (x min = 0.008), E665 [25] (x min = 8 x 10~ 4 ) and 
SLAC e~p experiment [26] (x m ; n = 0.063). Let us remind here that the lowest Bjorken x for 
the combined HERA analysis is by far the smallest: x m \ a = 5.52 x 10 -6 . This alone singles 
out HERA data as the best experimental sample to look for geometrical scaling. Since, 
basing on our experience with e~p data, we expect some possible systematic differences 
due to the projectile used in different experiments, we have analyzed - following the steps 
of Sect. 2 - each of these experiments separately. We have observed, indeed, rather large 
systematic differences between various experiments in the overlapping region of Bjorken 
x's (i.e. for large x) as far as (A mm ) and ((A m i n )) are concerned. However, due to large 
errors of (A m i n ) and ((A mm )) it is difficult to quantify the amount of violation of GS in 
these experiments. As far as E665 data is concerned, we could see qualitative agreement 
with the combined HERA data, although E665 data suffer from very large errors for (A mm ) 
and ((A m in)). More detailed analysis and comparison of DIS experiments [22-26] with the 
combined HERA data [17] will be presented elsewhere. 

It is worthwhile mentioning that our results for the best value of the exponent A are in 
a very good agreement with the quality factor analysis of Ref. [14] A = 0.321 ± 0.056 and 
with later analysis of the combined HERA data only [16] A = 0.31. In both cases, however, 
the analysis was constrained to x < 0.01 with some additional cuts on Q 2 and y which 
have not been applied in our study. 

Obviously our method can be applied to test different forms of scaling variable r 
which follow from the QCD non-linear evolution equations [14-16]. Finally, the question 
of the possible existence of geometrical scaling for the charm cross-section F2/Q 2 [15, 27] 
where the charm quark mass plays an important role may also be studied with the method 
proposed in this paper. 
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